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1.0 SUMMARY 

In this essay a detailed analytical and experimental study of 
anew, high speed, linear torque motor has been carried out. First, the 
need for such a motor is shown. This is followed by a general description 
of the actual motor. The equations of motion and impedance for the motor 
are developed theoretically and analytical expressions for all the constants 
appearing in these equations are derived. Experimental values of all the 
constants are then found and experimental runs are made of the phase and 
amplitude characteristics, the impedance, and the transient response to a 
step input. A comparison is made of the theoretical and experimental 
results including a discussion of linear and non-linear regions. The 
effect of the motor in the closed servo loop is then discussed as well as 
the effect of current and voltage feedback. Finally, the conclusions 


arrived at as a result of the study are stated. 


2.0 INTRODUCTION 
2.1 STATEMENT OF PROBLEM 

At present there exists considerable interest in and need for 
small size, light weight, high speed, linear, torque motors at the Applied 
Physics Laboratory, The Johns Hopkins University, Silver Spring, Maryland. 
These motors are used in closed loop, electrohydraulic servomechanisms 
(see Fig. 1). To properly analyze these servomechanisms, the operating 
characteristics of all components should be known thoroughly. It is the 
object of this thesis to conduct a detailed analytic and experimental 
study of one of these high speed torque motors. 

The motor was developed independently by the Applied Physics 
Laboratory and Curtiss-Wright for servosystem usage. A similar type 
was constructed by the Dynamics and Control Laboratory of the Massachu- 
setts Institute of Technology. Two types of motors were built for the 
Applied Physics Laboratory--one by Curtiss-Wright Corporation and the 
other by Bendix Corporation (Eclipse-Pioneer Division). No thorough 
combined analytical and experimental studies have been conducted up to 
present due to the lack of time. 

This analysis is conducted on the motor built by Eclipse- 
Pioneer since it is used in preference to the Curtiss-Wright motor in 
which the phase shift varies greatly with frequency because of large 
eddy current losses in the magnetic structure. 

2.2 GENERAL DESCRIPTION OF MOTOR 

The motor is essentially a polarized relay consisting of two 
magnetic yokes energized by a permanent magnet, a moving armature made 
of laminated iron, and two armature coils. The armature is pivoted at 
Point P, midway* between the laminated pole faces of the yoke, as shown 


in Fig. 2. One end of the armature extends out beyond the gaps and is 





WSINYHOSWOAYSS OIINVYGAHOYLO313 dOO1 G3SS010 V 40 WVYOVIC 40018 oS 


(39VLI0A v ATIVNSN) IVNOIS HONNS= °9-'@ = 


TVYNOIS LAdGLNO 
IWNSIS LADNI IWNYSLXS 


it tb 
aa w 


HOLYNLOV AINIVA YOLOW 
"2 YSAISNVYL INONOL 


YHAlal Id WV 


UJ 
D> 








C3dA1 AVISY G3ZINV10d) HOLOW 3NDYOL G33dS HOIH JO OILVWSHOS 2 Old 


xAI4 3YNLYWYEY = °P 


Ip 


Saicg 
AYNIVWYV NI ONIMOTS 
INAYYNOD WWILNSYS4SS510 


CSAANIVA YSSSNVYL) 
Qvo7l 


SYNLVWYY 





LANOSVW LNANVWYad 
JWOA DILINOVA 





connected to the load. 

The permanent magnet produces a constant flux, Do» which divides 
into two paths and generates equal flux densities at each of the four gaps 
when the armature is centered. If, under these conditions, currents are 
allowed to flow through the two armature coils then each current produces 
a flux. The difference of these two fluxes will be called the armature 
flux and will be considered to be produced by the difference of the two 
currents actually flowing or what can be called the differential current. 
This armature flux causes the armature to become polarized in such a way 
as to be attracted by one pair of diagonal poles and repelled by the 
other. There is a torsional spring restraining the armature motion which 
together with the design of the magnetic circuit causes the physical dis- 
placement of the armature to be proportional to the current in the coils. 
Reversing the direction of the differential current in the coils reverses 
the direction of motion. Torque is produced by simultaneously raising the 
current in one coil and lowering it in the other. The motor is operated 

with a quiescent current flowing in each coil. The motor operates with 
direct current power and is in principle similar to a Class A amplifier. 

An exploded view of the motor is shown in Fig. 3 and the flux 


paths for the armature centered and off-center in Figs. hA and LB. 
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3.0 THEORY 

In this section the theoretical equations of motion and impedance 
will be derived and theoretical values or expressions for all the constants 
appearing in them. 

3.1 DERIVATION OF EQUATIONS OF MOTION AND IMPEDANCE 

3.11 Equivalent Circuits and Equations of Motion and Impedance 

3.10. Linear relationship between force produced by the motor 

and the differential armature current. 

The first thing to be found is the relationship between the force 
exerted by the armature arm and the differential armature current. The 
flux in the armature, Ds is proportional to the difference of the currents, 
i, flowing in the armature coils, i.e. $= a i, where a is a constant as 
long as the iron is not saturated. In Fig. 2 it can be seen from the geome- 
try of the relay that the gap lengths gy + Bo always equal B34 E),° The 
permanent flux, oy emanating from the permanent magnet is considered con- 
stant. It will be assumed that as long as the movement of the armature 
is small and the permeance of the gaps low, the flux through the armature 
from the permanent magnet is negligible. 

The force at any gap is given by the following formula according 


Coumes. 1, p. 202. 


2 
F-= BS where F= force in lbs. 
«(on B= flux density, kilomaxwells/ sq. in. 
S= area of gap, sq. in. 
But, B= @/S where Z,is flux at any given gap 
so r= (Al's 
[S| 72 


If — O and the armature is centered, there will be no force 


present to move it away from canter, i.e. Pc Py and co B in Figses. 








If $,~ 0, then the total flux at each gap will be (see Fig. 2) 
At gy 2, = Bf2 + p,/2 
Bo = 8/2 $,/2 
g, 8, 8/2 - 2/2 
g), $= 0/2 + $,/2 
The force at each gap will then be: 


FL = (6/24 6,/2)° 8 
“i (ie 


FL, = (f/2- 9,/2)2 s 
2 San = 


hry 
Ww 
tl 


( B/2 - 8/2) § 
32 ie 
F, = (6/2+ B,/2)* § 
pee ee a 2 
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FIG. 5 FORCES ACTING ON ARMATURE, 
The forces acting on the armature are as indicated in Fig. 5 
(all attractive). 


Writing the torque equation for the armature, 


ae = a + ati 
FLX cea Fox ee F(X i) 


or Ft ee POT ae Ae) 





Substituting for Fy» Fos B3) and By the expressions just derived and at 
the same time putting in the numerical values of X and Y, namely X=0.781 


in.and Y = 0.425 in., one finds that 


rye Be 


64.78 
or FL= b 2. where b = g., 
Gliefs 
But g =al 


therefore, F,= abi-= Kyi where Ky = ab 
Ky will be known as the electromechanical coupling constant and will be 
discussed in more detail later on. ‘It is now seen that if the assumptions 
made in this derivation hold, the force is proportional to the differential 
current. 
3.112 MECHANICAL SYSTEM 

Now, considering the motor as a mechanical system, any force 
produced by the current in the armature coil must be numerically equal 
to the forces due to the inertia of the armature plus the moving mass 
of the piston and piston rod, to viscous friction, and to an effective 
spring force. By effective spring force is meant the combination of the 
force due to the torsional centering spring and the force due to permanent 
flux in the armature as a result of displacement. This mechanical system 


can be represented diagrammatically as follows: 
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FIG. 6. DIAGRAMMATIC REPRESENTATION OF MECHANICAL SYSTEM 





The symbolic diagram of Fig. 6 is: 
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FIG. 7. SYMBOLIC DIAGRAM OF MECHANICAL SYSTEM 


Writing the differential equation for this system, 


F—- Md°x, f dx, kx De 
dt dt M = 

f - 

k = 

xX = 


Since ie Kyi 





2 
Fy = Nei —- wes + £ dx |. kx 
ait dt 
3.113 ELECTRICAL SYSTEM 
The electrical system for the 
we Pig. o. 
a7: oa 
R= 
ir Ss 
| a 
eC ZL 3 Kp = 
| ie am 
| AA) 
pee heat Pe 
FIG 8. ELECTRICAL SYSTEM FOR THE MOTOR 


The output voltage of the last stage of 


sistance and inductance of the armature 


force produced, lbs. 


total effective moving mass, 

lb sec®/in. 

viscous friction, 1b sec/in. 

effective spring coefficient, 
lisyin. 

armature displacement, in. 


motor may be represented as 


resistance of coil, ohms. 
inductance of coil, henry. 
back eom.f., volt sec/in. 
voltage feeding armature coil, 
volts. 


the amplifier must feed the re- 


coil and also overcome a back 


*! 





ee 


electromotive force due to the velocity of the armature. An equivalent 


Circus Crawling. 9. 
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FIG. 9. EQUIVALENT ELECTRICAL CIRCUIT FOR THE MOTOR 


The voltage equation, by Kirchhoff's law, is: 


Keg = Ri; Ldi + Kpdx where R=rt Ry 
dt dt 


3.11) EQUATIONS OF MOTION OF THE MOTOR 
With these equations 


i Madox 4 fdx 4 kx 
dt ede 


and Hey = Hi + Ldi + Kydx © 
dt dt 
the output displacement per unit voltage input to the motor and also the 
impedance may be determined. Assuming that the system is initially at 
rest and using the Laplace transformation technique, one may write the 
following equations: 


(Ms° + fs +k)xX — Kuyt eit) 


(Kgs) + (Ls + R)E =e, 


A bar (-) over a letter indicates the Laplace transform of that quantity. 


Solving for X and i by means of determinants, 





0 hy 
ae R 
7 Me, Ls + Ky gp 
SL ae eee ee. Se 
Ms + fs +k “hy ieee fs +k) (Ls + R) + Kykps 
Kgs Ls + 





IMs? + (fL+ MR)s° + (kL + £R + Kpky)s + kR 








Ms° + fs +k 0 

ae Kps Le a (Ms°+ fst k eg 
IMs? + eu k = i IMs? + CaGs MR)s* + (kL + f£R + Kp ky )s + kR 
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The Laplace transform of the displacement per unit volt input is then 


x Ky 


— 


OR > fu ee 
pe IMs? + (fL+ MR)s- + (kL + fR+ KK s + kR 


= upon dividing 
s- + (f/M+ R/L)s* + (k/M + £R/LM + Kpky/LM) + kR/LM 


numerator and denominator by LM. 


The Laplace transform of the impedance, Z, is e/i where & =e, when the 


plate eee) 3 is zero, or in other words, one has & zero impedance 
source. 


Zo Pee ae | 53 + (£/M + R/L)s® + (k/M + £R/IM + Kpky,/LM)s + KR /in| 


z s° + (f£/M)s + kM 





where R now equals only ft 


The Laplace transform of the displacement per unit ampere is found by 


13. 





hie 


dividing x by i and this reduces to 


Ky 
(Ms¢ + fs + k) 





pet | > 
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3.115 MORE EXACT EQUATIONS OF MOTION 
It should be noted that the equations of motion derived in 
3.11) are simplifications of the true situation. No attention has been 
paid to the characteristic of the connection from the motor armature to 


the piston valve. A symbolic diagram taking this into account is shown 
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FIG. 10. SYMBOLIC DIAGRAM OF MOTOR TAKING INTO ACCOUNT ELASTANCE OF 
CONNECTION FROM MOTOR TO VALVE. 


spring coefficient of valve. 
spring coefficient of connection. 
displacement of piston rod. 
viscous friction of piston valve. 
mass of valve and rod. 

mass of armature. 
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Writing the equations of motion, 


peg = (Rot rp)it L di + Kp dx 
dt dt 


and Kyi = M,d°x + kx +Q 
at2 
Where Q=k, (x, — oe MX t sae a Ky 
from which x — : Mi + fa ie | 
K} 

A S*ngle dot over a letter indicates the first derivative with respect to 
time of that quantity and two dots the second derivative. 
Taking the Laplace transforns, 


HB, = (Ro+ ie Ls)i + Kpsx 


8 
= (Migs? + k)x + (Mje* + fys + kx) = (Mgs¢+ k — kX +k xy 


— 2 = = = 
ae (Ms + k k, )x + ky xX, 
By 
and X =_1 (Mjs@+ fys tk, — k))x, 
Ky 


Lerten = eet rye 
Then /le = (R+Ls)i+ K_sx 
/ g B 


or substituting for i and X 


{ ~ 
HE =(R+Ls) 1Ms*+ k —ky) J-1 Mys*+ fyst ky - ky) H+ kk 
Fe La | 
\ J 
Shes (ipa ec eaiee ek), 
Ky 
Dividing through by xy and then taking the reciprocal, 
ai 1 
nae ee. eRe ae nt) eee 
M&g (R+Ls) 1 geek ly (Ms F Aye ya} + -Kgs (Ms ¢fys + k,-ky) 
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Mi LM |e 





Multiplying the numerator and denominator by Kyk), we have 


K k 


= M 


SPST Tee 
keg ky + Ls as“ + k — ky) {-1 (Mys*+ fys+ ky — “e v 
k i 
io | 


— ones 


- Kk (M,s° + fys cake eo 
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a ~ By 
Aso+ Bett C's? + Dist t+ E's +F* 


where AX = Mam L 


&® 
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M fy R ate M tay 
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When k, is very large, i.e. k} ——> oo 
a= 0 
oe = 0 


C7 =-L(M +M ) 
Vv a 
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D R(M de M.) Lf 

py fh Jae) — Kaige 
By Sk k)R 


Substituting these values into the above equations and letting 


M = My+M, and fy=f, 


= Ky 
IG en Ba ee ile 
Ke, IMs-+ (fL + MR)s° + (kL + fR + Kak, )s + kR 





which is the same equation for displacement per unit input voltage derived 
in 3.114 for the simpler system assuming infinite elastance for the con- 
nection from the armature arm to the piston valve, an equivalent effective 
mass equal to the sum of the mass of the armature and valve, and that the 


viscous friction acts at the armature rather than in the valve itself, 


3.12 Derivations of Mathematical Expressions for Constants Appearing 
in above Equations and Calculation of Set of Specific Values for 
Motor under Study. 


In this section the various constants appearing in the equations 
of motion and impedance will be derived mathematically. Later these theo- 
retical constants will be compared to those found experimentally. 

3.121 Effective Mass M 

As has already been noted the mass of the armature must be 
accelerated. Since this motor is to be used as a part of a system and it 
is the object of this thesis to determine its operating characteristics in 
that system an adjustment must be made in this term to take the moving mass 
of the piston rod and valve piston into account. For the purpose of sim- 
plification the mass of the armature can be combined in an equivalent mass 
with the total moving mass of the valve piston and rod. 

i= Moo Me 

M.= mass equivalent of armature 


M,= total mass of valve piston and piston rod 
(measured and found to be 17 grams =0.0375 pounds) 


M = total effective moving mass converted to position 
of valve 


To find the equivalent mass of the armature, the moment of 
inertia of the armature was calculated from its dimensions and shape 
using the fact that the moment of inertia of a body is equal to the 


integral of the elemental mass multiplied by the square of the distance 








from this elemental mass to the axis about which the moment of inertia is 
sought and the parallel axis theorem. The calculation for I,, the moment 
of inertia of the armature,will be found in appendix 8.1. Theoretically, 


ee ene gm-cm?. Therefore the equivalent mass is 


a 18.1/(3.57)° = 3.78 gran or 0,00835 pound. 


re 
r = 3.57 cm is distance from pivot point P to point of 
application of valve force. 
Now the total effective moving mass, M, can be found. ee 


M = 3.78 + 17.0 = 20.8 pm or 1.188x1074 1b sec2/in. 
3.122 Effective Spring Coefficient k 
For static equilibrium the forces on the armature arm must be 
equal, i.e. FL = k 5x = Kyi + kpx 
where aro actual spring rate of torsional 
centering spring. 
Ke = anti-spring rate due to the permanent 
flux through the armature when armature 
is off center. 


Ky = electromechanical coupling constant 


The force, kpx, due to the permanent flux flowing through the armature 
when it is off center (see Fig. 4B) causes the armature to move farther 
off center. This force plus the force due to the flow of armature current 
must equal the force due to the centering spring for static equilibrium. 
The displacement sensitivity per ampere, x/i, is found as follows: 

kx = Kyl + kpx 

x/i ee (ko — ke) 
UK ee ke) is set equal to k, k may be looked upon as an effective spring 


coefficient. 








In order to find the value of k necessary to get mechanical 
resonance at, for example, 210 cycles per second, use is made of some well 
known relationships for second order systems. (These relationships may be 


found in any standard text on servomechanisms, for example, Ref. 2.) 


Oo, = fk w, = angular velocity at undamped 
M natural resonant frequency, 
rad/sec. 
k = effective spring coefficient, 
lbs/in. 
M = mass of system, lb sec“/in. 
a 
or k= Man = (1.188x1074) (1320)2 = 207 lbs/in. 


This particular resonent frequency was chosen because of the effective 
spring rate found in the experimental section and it must be borne in 
mind that it is only the true resonant frequency for a system having no 


damping. 


bele5s Viscous Friction f 
Assume J = 0.1 as a reasonable estimate as to the amount of 
damping present in the system, where i is defined as the ratio of the 


actual damping to the critical damping. 
=F = [ and 2)u = f/M 
M 
These relationships {= a ee f = viscous friction coefficient 


a oMon 2Mk 


or f= (2 Nk 
a = 
=(0.1)(2) .1.188x10 L007 S05 sae 





sec/in 
A useful parameter for calculations was found to be f/M which for | =oe 


was 


oe — 26h/sec 
M  1.188x1072 








TABLE I, VALUES OF £/M CORRESPONDING TO DIFFERENT VALUES OF f 


cme ee ee re ee 


* £/M 





| 0.018 7.5 
0.100 264 


0.113 300 | 
0.300 800 


The valve of Jy will be used as a parameter to show the variation 
of impedance and output displacement per volt as the amount of damping in 


the system is varied. 


3.12h Inductance L 


L— Nf, _ 7650 x 6560 __ 8. henries 





ix 108 0.060 x 108 
N = number of turns on armature coil 
a= armature flux in maxwells 
i = differential current in amperes 


The values of $, and i used to determine the theoretical value of inductance 


were taken from Fig.55 in section 8.2. 


3.125 Back Electromotive Force Ke 


This method of determining the back electromotive force (e.m.f.) 
is due to Mr. R. B. Higley and Mr. J. Wheeler of Curtiss-Wright Corporation 
(see Ref. 3, p. 17). The leakage flux, $,, however, has been determined in 


a different and more rigorous manner. 








x = valve displacement, in. 





y = armature displacement 
in alr gap, in. 


FIG. 11. RELATIONSHIP BETWEEN VALVE AND ARMATURE DISPLACEMENTS 


t+ can be seen from Fig. 1l that x W_ 


y 
1.06 0.781 


or oe OL = 0.6 
1.06 . 


When the armature is centered all gap lengths are equal or 
8) = 8) = 83 = &, = &- 
When the armature is not centered, g5= 3 = 8 — 0.6x 
g, = g, = Ct 0.0m 


The permeance at any gap with the armature centered is 


(P — 4S = permeability 
g So = area 
= gap length 


with the armature off center, = A, = 49 
“Pea. 6x 


ga 


0O.6x 





Let us draw the magnetic circuit for the flux of the permanent magnet 


(see Fig. 12). 


R,, Ro, R3, and), are 


reluctance of air gaps 
lS 2). ane 






ve 
PERMANENT 


y, is reluctance of 
Macner 


armature. 





FIG. 12. MAGNETIC CIRCUIT FOR PERMANENT MAGNET FLUX 


An expression for the flux, Drs flowing through the armature due 
to the permanent magnet can be found by solving the three magnetic loop 


equations. These equations are: 


Qyt+ Ro) hy — Ry go’ — Ro $3’= NI 
—Ry byt A + R3t+ Ky) f° = Ire ea, = 
ae Oy — K, fo’ + Qot+ Ry, + Ry) f3’ =0 
The flux we are seeking, Dr» is equal to ies P, 


Actually solving we find 


J, _ NI @8\-— @oR3) 
R, A3+ Ri) Rit R2)+- RM Ait Lot RiK0(K3+ Ai,) 


or, in terms of permeances, remembering ic Fi and iss F 


%, — NI Cz i: +n) 
afar a) ale* a) + aale * a) 


LL 





Dividing numerator and denominator by i + 4) 
/ 3 


bh = We -z 
e(e ta)? 


Putting the fractions over a common denominator and multiplying numerator 





and denominator by ee P 3 


o, — NIA, - &) 
E+ 
ge +e 


L. 





But we may justifiably assume, since the reluctance of the path through the 
armature is approximately zero, that the permeance, , is approximately 


infinite or certainly very much larger thant, + F 


Therefore, P+ F, <0 


Pr, 
as ae 


Substituting for F, and ta one has 


ee NI 14S 2(0.6)x | 
g°-0. 36x? 


CURA Seem x 


This checks with the expression given by Higley and Wheeler in Ref. 3, but 


solving in this manner requires only one assumption, namely, that EYEE. 





It should be noted that one can reduce Bo which is desired small, by either 
keeping x very small or by increasing g, the gap length, with the armature 
centered. 

To get the back e.m.f. due to oe use is made of the fact that 
2, is a function of displacement x, and x is a function of time. Therefore, 


Dy is a function of time. 


E =~ Nagyic~® volts where N is number of turns of coil 
at g, is flux in maxwells 
E = —Ndgy dx 1075 
ax Jat 
B= K, dx where K,=—N ddr, Mere volt sec/inen 
dt dx 


If we differentiate ¢, = NIKS i<2 x with respect to x, we get 
é Be 3676 


= = 
df, — NIuS 1.2 | (g@-.36x2)-2 4 2(.36)x2(g2~.36x2 -2| 
= sh g x doy 2( 30) x“ (eo—, Ox) 


_NT ys 1.2 | (g*+ .36x°) 
2 (g?-.36x2)¢ 


Substituting this into our expression for Kp» we get 


ee (NaS pec gt .36x*n0-8 volt sec/inch 


B 
2 rc ~.36x?)2 





To use this formula the NI of the permanent magnet must be determined. 
The method used to determine this was as follows. The flux flowing in 
the armature due to the permanent magnet with the armature centered was 
assumed negligible. Also there was assumed to be no leakage flux. With 
these conditions plus the further assumption that the reluctance of the 
iron path is negligible compared to the air gaps, all the magnetic poten- 


tial difference of the permanent magnet will appear across air gaps 1 and 2 


(see Fig. 2). 


~~ 
7 = R Bp f == magnetic potential difference, 
gilberts 
A == reluctance, gilbert/maxwell 
(Rp eg ~p = permanent flux, maxwells 
a8 A. = permeability, maxwell/NI-cm 
By = flux density in any gap due to 
2o= BYS permanent magnet, gauss 


Substituting these expressions for R and Do into iS 


raft = 


But - gilberts = 0.l77NI ampere-turns 





SO, NI — 29B, = 2(0.145)h170 —766 ampere-turns g = 0.115 cm. 
Sie nee (measured ) 
Ob 0.bur{1. 259) Ato = 1.259 maxwell/NI-cm. 
B, = h170 gauss 


(measured ) 


By = 4170 gauss was determined by use of a General Electric fluxmeter 

Model 32C2h8 which was calibrated at the National Bureau of Standards by the 
author for this purpose with the assistance of Mr. I. L. Cooter and A. R. 
Lindberg of the staff there on 2 August, 1951. Substituting in the equation 
for Kp and letting | 


aa Kpmin = 9.01 volt sec/inch 


x =0.032, Kpmax =12.6 volt sec/inch 


3.126 Electromechanical Coupling Constant Ky 


Ky, the electromechanical coupling constant, may be defined as 


follows: 
Ky — F (force in lbs) 
i (differential current in amps) 
But F (~p+ By)? from section 3.111 upon making the proper 


2888 


substitutions in the formula for force. 








2 
(Dot Ba) 
Therefore, Ky= 2008S where 2%, = 13.9 kilomaxwells (measured) 
41 
$, = 6.56 kilomaxwells (Fig.55 ) 


S = gap area, 0.262 sq. in. 
i = 60 milliamperes 


Upon substituting the numerics and solving, Ky is found to be: 


KM — poe ens — 92,8 1bs/amp 
{2BBY (0.262) (0.060) 
3.2 PLOTS OF THEORETICAL IMPEDANCE AND PHASE AND AMPLITUDE CURVES 
3.21 Impedance Curves 
In section 3.11 it was shown that the Laplace transform of the 
impedance was: 


2 L[s3 + (£/M + R/L)s2 + (k/M + £R/LM + Kpky/LM)s + kR/LM] 
. a4 fs/M + k/M 


s 
or multiplying through by L, 


Z_ us + (£L/M + R)s? + (KL/M + £R/M + KBKM/M)s + kR/M 
3° + fs/M +k/M 


Making use of the properties of the Laplace transform, since s=vVv +jw 
with7=0, S = jw, Now,writing the impedance as a function of frequency 
by substituting jo for s, one has: 


2( jw) Ae ( 5a) ss L(4w)? + (£L/M + R)(40)* + (KL/M + £R/M + KpKy/M) jotkR/M 
(30)© + £(j0)/M + k/M 


orseparating the real and imaginary parts 


2( jo) [R(k/M-w7) 10°] 4 joo [L(k/M-w) + £R/M + Kk /] 
(k/M—e9") + §(£/M)oo 


Rin this case equals R, only. 


From this equation tor impedance plots have been made in Fig. 13 of Z as a 


function of frequency. Also the variation of Z as the viscous damping term 
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TABLE IIT CALCULATION Dala FOR THBCRSTICaAL InreDacCi CURVES 


49200 


L/ a = _ £7 =47.5 £/.d = 300 
We kee re 10 OS re 0 
f Zo It, 2 © ei g a Q 
cpS Yrad/sec’ ohas degrees : ohas de,rees ones degrees - 
| 
0 0 | SOD Jao 1800 Oe?) hago Oo 
eX) 6 Pe e970 52.0 S029 Doon 
40) aor. ola 69.4 
60 376 7450 76.0 Lone deel: 
30 Soom 89271 4 7956 
LOO feo, Laeldo 81.2 Pao) oy See Le 7 Co 
120 794 (oO mee. O 
149 os.) 1y400 34.7 19490 fie: 
1690 Loos 22d) 895.7 
3) io) wooo)? S640 a0) Boe eat iia 
HE) Mest Me Ge650 (38756 e730 6329 
Bo) to awos6le fs. or. 638800 esc 44909 ou 
<1LO Mec te of =9).79 148000 a SEO Goa a 
ae ee SoZ 0a 32S CMC = 55.0 GOGO oe 
eo) 144% 75990 -76.8 Hoe. =13.6 Leno! Le 
W~40 eae 7 0649 -650.5 eee +61, 7 26109 AD eo 
w44 nee Bice Done. toe 
Zo) moto « Eseo0 Wo. 
BO) LOGO eo Oe eu. 3 16500 69.4 
aU) 1760 (oOo eae 
BIO £5350 ole he) ae elo 34.0 22100 BOL 
540 z140 ke Foo Sore 
560 ~2£60 °° 69400 386.5 
49 oS POSS see Secu 26.) 342009 36.0 
ae EOe Oo = ee 





is varied. Calculation data for these theoretical curves may be found in 


Table II. 


3.22 Frequency Response Curves 
In section 3.11) it was shown that the Laplace transform of the 
displacement per unit volt was: 
—% _ Ky/LM 
ME st + (£/M+R/L)S°+ (KAP + ER/IM+ Kpky/IM)s+kR/IM 
or writing this as a function of frequency as explained in 3.21 
x( Joo) Ky/ LM 
Fee — (kR/IM - (£/M + R/L)o*|+50[k/M + £R/LM + Kyky/LM - 02] 
Plots have been made of XY/Le p versus frequency for various values 
of |. In Fig. 1h) curves are plotted considering the tube as a zero impedance 
source, i.e. fo The extremely large phase lags prohibit the use of 
such a source for driving the motor in the servo loop outlined in Fig. l. 
In Fig. 15 curves are plotted considering the plate resistance, Ty? to be 
equal to 60,000 ohms and it is obvious from a comparison of the change in 
the phase and amplitude characteristics why pentodes with their high plate 
resistances are used to crive the motor. Calculation data for these curves 
may be found in Table III. 
Plots have also been made of the displacement per unit ampere, 
x/i, versus frequency assuming a constant current source in Fig. 16. 


Calculation data may be found in Table IV. 


3.23 Transient Response 
To secure the transient response characteristic of the motor a 
step input voltage, K; is applied. That is, 
a f(t) =(0, t<O 


ee 10 
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TABLE WV CALCULATION DATA FOR PLOT OF OUTPUT DIS LACEMENT 
PER AWPERE VERSUS FREQUENCY 


a ee 


“yoo Fo 


eee ae | i 
0 TAT | 5 2 ater O ; 
20 eso 6 Poa, Oe 
MGM 0.567 1 =1.5 4. 3.553 in 
moet O.%2sy) 2.5 9) O.7ea |i 16.41 
| 189 | 1.65 | Cau ie alan 36.4 | 
fee | ees |) xo.5 91 1285 67.0 | 
plooey| ews | 90.9 | ade 99.0 | 
PO 8 4.15 + 155.3 § 75 Plo | 
moo) 2.19 | léd.6 ies Mo le8.e 
“40 i 1.39 | 174.5 1.08 140.1 | 
po.) # 92366 | ys.1 pede og 151.7 
eoO “Portis | i77-e GA. 400-7. 16.7 . 
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Taking the Laplace transforn, ue, = K/s 


Substituting this for “ue, in section 3.l1h, we have 


eo Km K/s 
IMs? + (fL+ MR)s¢ + (kL +f£R + Kpky)s + kR 


| 


| 


or 


Ky J 
5 [ IMs? + (fL +MR)s* + (kL + £R oe 


Dividing the numerator and denominator by kR, we have: 
e Ky/kR (K 
S [ (N/R) s? + (£L/KR + M/k)s* + (L/R+£/k + KpKy/KR) s +1 i 


| 


and substituting in the numerical values (values listed in appendix — 


Table XVI under motor B) 
x eeorc alos Ke 
ae = = = 
3(2.29x1072 953 + 76.75x10~s2 + 6.07x10~s + 1) 


Factoring the cubic in the denominator, 


Eo 9.64 x 107° ie 
saa oon) a - i .. eS 
s| s+ 2.95x10° || (s+203) + j1190]| (s+203) - $1190} 


If now the inverse Laplace transform is taken, 


-2950t -203t 


x= 9.61:x0"° [1 .000 - 0.162e +1.000e sin(1190t 123.1") | K* 
At t—0O, x —9.64x107° (1.000 270.162 = 0.830)K 
= 
and at t=0cco, x=9.6x107>K* 
If K=1l, that is a unit step voltage is applied, 
x = 9.6x1075-1. 56x1072990F +. 9. 64x107e7 203 tgin(1190t-123.1°) 


A plot of the transient response to a step input of one volt is shown in 


Fig. 17 and the calculation data in Table V. 
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4.0 EXPERIMENT 

In this section the actual experimental techniques used to 
determine the constants derived analytically in section 3.0 on theory 
will be described. It will be noted that two motors, called A and B, 
were used. When the analytical section was begun it had been planned 
to use Motor A throughout both the theoretical and the experimental 
phases but in the course of the determination of the phase and amplitude 
characteristics the laminations composing the armature of Motor A separated 
and therefore a different motor had to be substituted. This new motor has 
been called Motor B and differs from Motor A only in that it is a separate 
physical entity subject to the usual variations in manufacturing processes. 
For the comparison of theoretical and experimental emits a new set of 
theoretical curves for Motor B was computed. The values of Motor A have 


been left in the report for the purposes of comparison and to retain the 


analytical studies already completed. 
l.1 EXPERIMENTAL DETERMINATION OF ALL CONSTANTS 


lll Effective Mass 
The moment of inertia of a body can be determined experimentally 


from the following formula (see Ref. lh, p. 72). 


dee Tew MaG where I, —moment of inertia of 
<n armature in gm~cm?. 
a T = period of oscillation sec/cycle. 
ow == suspension width in cm. 
M, —mass of armature in grams. 
1 = =length of suspension in cm. 


G =acceleration of gravity, cm/sec¢ 





LMM A 


Li 


ARMATURE 


FIG. 18. BIFILAR SUSPENSION OF ARMATURE FOR EFFECTIVE MASS DETERMINATION 
The’ armature was suspended as shown in Fig. 18. It was then set 
into motion with a very small amplitude of oscillation about its vertical 
axis. The period of this small oscillation was determined with the aid of 
an electric timer. The experimental values to be used in the formula for 


inertia were determined to be as follows: 


1 
T= 20 cycles =0.755 sec/cycle 
15.10sec 
wz 2.10 cm. 


M,a= 3.7 grams 
G = 980 cm/sec? 
1 = 56. cm. 
Substituting these values into the formula, 


I_—0.755° x 2.10% x 3.7 x 980 —h8.h gm-cm*, 
oe es 
ly 56.4 


x 
in order to determine the effective mass of the armature at the point of 


application of the valve force use was made of the fact that the moment 


a77e 





of inertia is equal to the product of the mass and the square of the distance 
to the axis about which it acts. Therefore, 


M = h8.u = 3.80 grams 


e= lg 
“re (3.57) 





where r = 3.57 cm is the distance from the point of application 
of the valve force to the pivot point P of the armature. 


Me= effective mass of armature. 


The total effective mass of the system,M, as defined in section 3.121 is the 


sum of M,e+ My 


or M = 3.8 +17.0= 20.8 erame (1.188 x 10741b @ec2/in) 


4.12 Effective Spring Coefficient k 
The method used to determine the effective spring coefficient, k, 


was as follows. A quiescent current of 2h milliamperes was allowed to flow 
in each of the armature coils and then standard weights were added to the 
arm of the armature at the point of application of the valve force. The 
weights were added very carefully and then removed in inverse order of put- 
ting them on to preserve the hysteresis loop. A Brown and Sharpe dial in- 
dicator was mounted on the motor in order to read the inches displacement of 
the armature arm. This dial indicator read to the one one thousandth of an 
inch from zero up to thirty thousandths. After the plotting of the experi- 
mental points, a straight line was drawn through the middle of the result- 
ing hysteresis loop and the slope of this line was determined to be 207 
lbs/in for Motor A. This was taken as the anne of the effective spring 
coefficient, k, (see Fig. 19). 

The value of the natural (i.e. undamped) resonant frequency may 


now be determined. | 


—_Wk—|| 207 —1320 rad/sec. 
ae it |e raerost 
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or the frequency is f = 1320/277 =210 cycles/sec for Motor A. 

Table VI contains the experimental data for the effective 
spring coefficient for both Motors, A and B. The data for Motor B 
are plotted in Fig. 20 and the slope has been determined for dis- 
placements of the armature to either side of center. For Motor B, 


k = 168 los/in. and f = 190 cycles/second. 


4.13 Viscous Friction f 
The value of the viscous friction term, f, was determined in 
the following manner. An impulse was applied to the armature and a 
photographic record was taken by means of a Miller recorder. The value 


of y was obtained from the record as follows. 


tim e 





FIG. 21. SKETCH OF TYPICAL TRACE ON MILLER RECORDER 
If in Fig. 21 we allow Q to be the logarithmic decrement of the 


envelope of the decaying oscillation, then 


= 1n y/y2 where y and yp are two consecutive values 
Z\ t of amplitude of same sign read from 


record. 


At is the time interval separating 
these amplitudes. 


It ig well known that (= [Wy 


therefore, { — ic where 6),=undamped natural resonant frequency 
() Vis : 


n 
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EXPERIMENTAL DATA FOR EFFECTIVE SPRING COE?!FICIENT 


- mete 


MOTOR A 


ee i ee . 


: 


{ 
b 


Left of Center - 
INCHES 
J sO007 4 
POOL] 0019 
265 Oe 
Jeo 5 BOO 
ON Pe) Ls) 
wee onal 
~woOrSsS ~OLIL 
MOZOS Oe. 
0 (aya. Aye FaKOYs 
Joy 3243 
TOO u wee 
,0256 259 | 
mozoL s0E6S § 
eZ z 


HO 


Uzoa8 


-_——_ 


ee 


0 
JO0I 
O79 


an) JeKs 
«JLES 
sO SC 
nO 
0240 


° 


»0Z99 


0269 
Jet? 


_—— --aoe we mee ee 


Left of Center 


ES OP ETE Sete ON ee a 





A rt a rt re et +) ee een “eens a eee neet ee. 


MOTOR 5 
KRigat of Center 
INCHES 
0004 0 J0ge 
FOOLS POO? ar 
sooSa7 0079 melerrs 
SOG POURS «J Les 
Bol73 . one 0189 
ALC, so205 Jee 
Ua426 WOZ oe «Jace 
iveoe ~J200 OL6F 
SOO nm Poors om te 
Foo 
Oats 





€# YOLOWN LDNEYNHOVIdSIG GHOLYRYY SNSYHA FAOYOT FYATIVA Os OIA 


27e" 


tpt 
an ie 
tate tt 


leat 
5 


ees. eieke 


eee ee 
tata. b 44 
-_ - + 


= 


pe dis 


| al: 


EE oe eee ee 


—- 








/ a where &,=damped resonant frequency read on ~ 
oscillographic record. 
(see Ref. 2, p. h9) 





oe TOR == =e 
aaa 
Substituting, 
ee /n a 
GAG 


{\7—f= 
Multiplying numerator and BeLerdn ator by We 


Tie: In $e VT =I™ 


Qa os & 
Squaring and separating terms, , 2 
2 ie /n % ; y° 
Vv — — eo aes 2, AE 


Rearranging terms, 


/n > e 2 
ig (Hi ie (/n a) ) 
_ L2 e) (a, At) + Un ¥ ae 


/n #— 


TG Get) == (in $-)? (In $y 


By means of this formula the values of J were computed from values 


Finally, 


of a2, AT, yy, and yo given by oscillographic records. A plot of the values 
of amplitude on semi-log paper should have given a straight line, but this 
did not prove to be the case, showing that there are non-linearities in the 
system and that the damping present is not pure viscous damping. Average 
values of all the runs made (about three per set of given conditions) were 
used to calculate the final value given. It should be emphasized that these 
values are subject to some error due to the nonlinearity mentioned above 

and the difficulty in reading the record. Table VII shows the values of 

q for various conditions. Runs were made with and without the amplifier 
connected to the motor and with and without the 0.01 microfarad condensers 


across the motor coils. 


=a 
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TABLE VII. VALUES OF CRITICAL DAMPING RATIO, 






CONDITIONS 






No amplifier, no condensers 






No amplifier, with condensers 










Amplifier, no condensers 0.13h 
Amplifier, with condensers 0.257 
Voltage across coils, with 

condensers, no pick-up Oeleo 










Voltage across coils, with 


condensers, with pick-up 0.229 








Knowing J , it is an easy matter to determine the viscous friction 
of f, since f =2JuW,M. Taking a value of) =O-1, for the case of no amplifier 
and no condenser, 


f=—2 x 0.1 x (190 x27) x 1.188 x 10-%=0.0283 1b sec/in. 


This is the value of viscous friction for Motor B. 


4.1 Coil Data--Resistance and Inductance Rg and L 


The resistance, R,, of the coil was determined by means of a 
Wheatstone bridge to be 1807 ohms for one coil and 179) ohms for the other. 
A value of 1800 ohms was used in the calculations for Motor A. For Motor B, 
the value used was 1750 ohms. 


The inductance, L, was determined as follows: 


Vz, = 140 volts, 00 cycle. 


V3 = 123 volts. 





FIG. 22. CIRCUIT USED TO MEASURE INDUCTANCE OF ARMATURE COIL 





The coil was placed in series with a resistor of 12,000 ohms 
resistance and a source of 00 cycle per second power as shown in Fig. 22. 
Voltages Vj, Vo, and V3 were measured. The current, I, flowing in the 
circuit can be determined since we know the voltage across the 12,000 


ohm resistor. 


I = 36/12K =3 x 107 amps. 


The magnitude of the impedance of the coil, Z., is the voltage across the 
coil divided by the current flowing through it. 

Var 2,1 = 123 volts 

Ze = 123/3 x 10°? = ki x 10° ohms 
The source impedance, Z,, is the voltage, V,, divided by the current, or 
with V, =1h0 volts 


Z. = 1h0/3 x 107? = 6.7 x 10° ohms. 


By drawing a vector diagram and actually plotting these values as in Fig. 23, 


one finds that 


X, I =115.0 volts 
Xy = 115/3 x 10-3 = 38.3 x 103 ohms 


in ==" 3653s 103/2r400 = 15.3 henries 





Va=AL 
FIG. 23. VECTOR DIAGRAM OF VOLTAGES 

The same procedure for finding inductance, L, was followed at 
60 cycles per second, a 60 c.p.s. source being used in place of the 00 c.p.s. 
source used above. The following voltage readings were obtained: 
V, = h5 volts 


Vo = 36 volts 
V3 =17 volts L =1h.1 henry. 





The value of inductance used for Motor A was 15.0 henry. Similar 
measurements on Motor B indicated that the inductance L was approximately 


20.0 henries. 


4.15 Back Electromotive Force Kp 


To determine the back electromotive force experimentally a small 
variable speed motor was connected to the torque motor through an eccentric 
cam. When the motor made one revolution, the cam caused the torque motor 
armature to move a certain distance x in simple harmonic motion. In other 
words, x = Awsinewt where A = maximum displacement of armature 

from neutral. 
& =angular velocity of motor 
The derivative of x with respect to time, or the linear velocity of the 
armature, is dx/dt = Awcosat. The voltage produced by the motion of the 


armature in the magnetic field is proportional to the velocity of the 


armature. 
i.e. V,=Kp dx/dt where V, = voltage generated by the 
& & 
motion of the armature 
or Kp =V at in the magnetic field. 
& Kp factor of proportionality 


The velocity of the armature in simple harmonic motion will be maximum when 


x 1s zero. Therefore, the voltage, Ves may be written as Ve = Vemaxcosat. 
Substituting this into the expression for Kp and at the same time substi- 
tuting for dx/dt, we may write 

Kp — V emaxc oset _ V max 





A @& coset AGW 
We can now find Kp experimentally by measuring the maximum or peak voltage 
generated, the maximum displacement of the armature from neutral, md the 
angular velocity of the motor. The voltage was measured by placing a 


Hewlitt-Packard voltmeter across one of the coils; the maximum displacement 


of the armature was determined by the eccentric cam; and the angular velocity 
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of the motor was found by the use of a stroboscope. 

Runs were made with eccentric throws of 0.00725 in., 0.011 in., 
and 0.0285 in. In each case graphs were plotted of the voltage generated 
versus armature linear velocity. These graphs are shown in Figs. 2h, 25, 
and 26. In order to determine the value of Kp from these figures, the 
actual slope in volts per inch per second must be multiplied by 1.1 as 
the ordinate is the voltage read on a high input impedance voltmeter or, 
in other words, the root mean square value of the voltage whereas the 
peak voltage is desired. The experimental data are tabulated in Table VIII. 
In Fig. 27 a graph was plotted showing the variation of Kp with the ampli- 
tude of the armature motion and on the same figure the theoretical values 
of Kp from section 3.125 as well. It can be seen that for Motor A the 
theoretical and the experimental results agree very well. Only the theoret- 
ical values for Motor B have been plotted. It was felt that since the 
theoretical and experimental curves for Motor A agreed so well and time was 
at a premium that a satisfactory value of Kp for Motor B could be read off 


the graph as long as the maximum amplitude was kept below 0.015 inches. 


4.16 Electromechanical Coupling Constant Ky 


To determine the electromechanical coupling constant, K,, 
experimentally, standard weights were hung on the armature arm and then 
the differential current necessary to bring the armature back to its 
center position was recorded. This procedure was repeated in taking off 
the weights and then the whole process was repeated for displacements on 
the other side of the center position by turning the motor over. This 
gave values of Ky for both positive and negative differential currents. 
To maintain the center position of the armature, a Brown and Sharpe dial 


indicator was mounted on the motor with the indicator set at zero. This 
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BACK ELECTROMOTIVE FORCE VERSUS ARMATURE DISPLACEMENT 


FIG «7 





gave the needed reference to enable the armature to be brought back to its 
center position each time a weight was added or subtracted to the armature 
arm. These data,as plotted in Fig. 28, gave a line whose average slope 

was determined to be 90.1 lbs/amp. for Motor A. The data for Motor B are 
plotted in Fig. 29, and Ky, = 100 lbs/amp. The calculation data may be 
found in Tables IX and X for Motors A and B respectively. Currents i, and 
i, are the currents flowing in each of the motor coils and i is the differ- 


ential current. 


4.2 EXPERIMENTAL DETERMINATION OF IMPEDANCE 
In order to determine the impedance of the motor experimentally 

a vacuum tube voltmeter was placed across one of the coils and a 27 ohm 
resistor was placed in series with it. By placing a second vacuum tube 
voltmeter across this known resistance, the current flowing in the coil was 
easily determined. The magnitude of the impedance was found by dividing 
the voltage across the coil by the current flowing in it. A plot of this 
experimental impedance as a function of alternating current frequency may 
be found in the section on the comparison of theoretical and experimental 


results in Fig. hh. 


4.3 EXPERIMENTAL DETERMINATION OF PHASE AND AMPLITUDE CURVES 

In determining experimentally the phase and amplitude charac- 
teristics as functions of frequency, the first problem encountered was 
that of measuring phase as accurately as possible. To do this a phase 


shift network was used as follows. 
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FIG. 30. SCHEMATIC OF PHASE FIG. 31. VECTOR DIAGRAM FOR 
SHIFT NETWORK PHASE SHIFT NETWORK 


A voltage, ej, was put in across terminals D' and F' and a voltage, Bo» 
was taken off the terminals A' and B' as shown in Fig. 30. The phase 
angle, %, is equal to 2 tan~1X, as can be seen from the vector diagram in 
Fig. 31. This phase shift aah es 1s discussed in Ref. 5, p. 30. Two 
matched resistances of 15,000 ohms each were used as R; and Rj'. A standard 
decade capacitor (Cornell-Dubilier Decade Capacitor, Model CDB-5) was used 
for the capacitance. 

Fig. 32 shows how this phase shift network was used to measure 
the phase shift between the input voltage to the amplifier and the output 
of the torque motor. The input voltage, e;, was fed into the amplifier 
and also into the phase shift network as shown in Fig. 32A. The signal 
from the amplifier actuates the motor and by means of the capacitor pickoff 
device a voltage representing the position of the armature arm is obtained. 
This voltage from the pickoff was fed into one pair of plates of a cathode 
ray oscilloscope. The signal, e,, from the phase shift network was fed into 
the other pair of plates. As shown in Fig. 32B, when the two signals are 
out of phase the trace on the oscilloscope will appear as a loop. When by 
means of changing the resistance, R,, in the phase shift network the two 


Signals are brought into phase, the loop closes to a straight line. Thus by 
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using the oscilloscope as a null detector in this manner, one is able to 
measure the phase difference between input voltage and the displacement 

of the armature by recording the value of R, necessary to close the loop 
at various discrete frequencies. The phase angle, %, is as mentioned 


earlier 2 tan“1x,. A photograph of the actual experimental setup used 


to measure phase shift characteristics is shown in Fig. 33. 

The capacitor pickoff deserves especial mention as it is a 
very sensitive device for changing the position of the armature arm into 
a voltage so that one may measure accurately its position. The circuit 
for the armature pickoff is given schematically in Fig. 3h and a photo- 
graph of the device as mounted on the motor is shown in Fig. 35. The 
pickoff is shown connected to the armature arm by a wire as it was con- 
nected during some early tests. However, in all the runs actually 
reported on the capacitor pickoff was mounted directly on the armature 
arm. The fan shaped piece moving between the stationary plates in accord- 
ance with the motion of the armature changes the capacitance, thus giving 
rise to the signal voltage. A plot illustrative of the linearity of this 
pickoff is given in Fig. 36. 

The amplitude of the armature motion was obtained from the photo- 
graphic record taken with a Model Hi, Multi-channel Oscillograph manufactured 
by the William Miller Corporation, Pasadena, California. An input signal 
was fed in by means of the signal generator to the amplifier and the motion 
of the armature of the motor was recorded on the oscillograph. Care was 
taken that the amplitude of the armature displacement remained within the 
linear range of our Kp factor. Also the voltage across the tube itself was 
kept above 65 volts at all times. The reason for this will be discussed in 
section 5.21. The Miller record gave a picture of the movement of the 


armature in response to the sinusoidal input. The experimental setup used 
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to measure the amplitude characteristic is shown in Fig. 37 and a sample ; 
portion of a record is enclosed as Fig. 38. 

‘The amplifier used to drive the motor is shown schematically 
in Fig.939. 

In Fig. hO the variation of phase lag versus frequency for three 
different initial values of amplitude is shown. It had been noted in per- 
forming tests on phase lag characteristics that the shape of the curve was 
easily reproducible but that the initial value of phase lag varied. Initial 
values of the amplitude of the input of 0.0025 in., 0.005 in., and 0.010 in. 
were taken and the resulting phase characteristics plotted. At low frequen- 
cles there is less phase lag for the larger values of sinusoidal input. This 
is due to some non-linear effect; perhaps Esterase which would cause more 
phase lag percentagewise for small displacements than large ones. As the 
frequency increases all the curves come closer together until beyond resonance 
the curves of the larger values of initial amplitude exhibit more phase lag 
than the smaller ones. A valid reason for the reversal at high frequencies 
is not known at present. 

The experimental values of phase lag of Motor B, both with and 
without the 0.1 microfarad condensers are shown in Fig. 1. It is seen that 
the addition of the condensers introduces an additional phase lag above 
resonance. 

In Fig. 42, a plot of the amplitude ratio (output to input) 
versus frequency has been made for Motor B, both with and without the 0.1 
microfarad condensers. In the same figure the theoretical values of ampli- 
tude ratio for the motor have been plotted for J =0Q.1. The experimental 
data for these curves may be found in Tables XI and XII and the theoretical 


value of amplitude ratio in Table XIV. 
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FIG.38 SAMPLE PORTION OF MILLER OSCILLOGRAPH RECORD OF 
ARMATURE RESPONSE TO A SINUSOIDAL INPUT VOLTAGE 
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EXPERIMENTsL DATA FOR PHASE 








TABEE 1 AND AwPLITUDE FR 
ee. ve B FREQUENCY 
CilAR AC : — ea aoc aes 
— a OF MOTOR B WITHOUT COADCNSERS ACROSe 
rr er ne 
| | a | Ko (2 tan “as pee | PATIO “ en xlO | 
“eps | ohms | ohms ‘de ae ty no 
ee} “ohms [de,rees” inches infvolt | 
Me | 2403 |31200 , "8.8" 0.41 | 220) eres 
18 | 1700 £0800 9.40.41 | 1.00 | 9.38 | 
29 | 1409 j15600 | lo.2 9.42 | 1.08 | 9.87 © 
40 | 1090 | 7800 , 14.6 0.44 | 1.06 | 12.0 | 
69 | 1000 | 5209 : £1.38 0.46 Pate | 10.5 | 
ao | 940 | 3900 , £7.2 0.50) lee | 1.4 
109 | 1090] 3120 35.6 0.55; 1.04 | le.6 | 
Veo); 1070 | 2600 | 44.8 0.62; 1.51 | 4, | 
Beem i220 || eso. 57.4 «On 76 me.) as 
; 163°) 1550 | 1950 | 76.8 90.37 we | ce. | 
130 : Bee lrcl A dese Gee erat «| e8.e | 
i99 | 6000 | 1641 ' 149.4 1.x4! 3.0. | 28.3 | 
Pace) .y fio00 (ess el. eee 2.c6 a | 
| 210 190990 | 1436 180.9 9.90 le.19 | 20.5 | 
| «20 | 100) 14ly 183.0 0.67 / 1.63 | 15.3 | 
ese ceo | 1500 4 19k amo, 44\8 1.97 i eo 
eco @ © 28) | 1200) 2062 0.32 | 2.777 730 | 
Meco F 339 | 1113, 210) 0.04; 0.584 | 543 | 
| 209 , 2301 1040 415.6 0.69 | 0.485 aac 
| 240 | 230} 919 419.6 | 0.13 | 0.216 Lg 96 | 
© 860 | 50 | B66 end. OVD | en “ae 
8 le ee 1.83 | 
| 450 | 380) 694 257.4 0.06 9.146 ae 
ee ee ee ee | 








TAbLE XII EXPERIMENTAL DaTA FOR PAASE AND AMPLITUDE FHE.U wwCY 
CiisnRACTEHISTICS OF MOTOR B WITH CONDEWSERS aCROS. 
wWOTOR COILS 


| po Bk LaweLstu¥k 
PF GR | Ae ten”) | nisee) | RATIO 
a aT Te 
ete | 500 im aS.8 | cc oe co 
[eee | 1600 250) |) ee Ogee 1 Lae 
47.5 | 1490 Boo | 2506 | Geeeo |) dots 
66.2 | 1400 | 4800! 32.6 | 0.385 | 1.28 
pom | igddm, O40 : 63.4 )"o440 | AES> 
104.9 | 1309 : $060. Gre | mcee! | eines 
lex. | 1400 | £690 | 57.6 | 9,610 | 1.33 | 
141.9 | 1630 | 2260 | 73.6 ; 0.799 | 4.4% 
159.0 , 2239 | £990] 96.e ; 1.019 | 3.11 
Meo | goo | 70)! Ieem I l.eo | 3.4% 
133.9 | 4400 ; 1740 | 136.8 | 1.a70 |! 3.e9 
132.0 | 6200 | 1669 | 159.9 | 1.020 {| z.14 
e19.9 | 49909 151) | 189.0 ; 0.855 | 2.65 
259.0 7 130 |. 1330 Peele | 0.705 i a: 
moeeo | 890d) | TaeOM zac.2 | 0.600, 1.35 
275.0 | 590 | 1160, 226.6 | 0.455 ; 1.49 
Po0l.3 | 670 | 10609 8244.3 | 0.245 | 1.6 
233.9 ; 3929 | 941 © 272.6 | 9.240 , 9.753 
oreo) | 1133." “ene ’ 28e.¢ 1 3.160 . 0.495 
AV. seo ee’ 300.6 | O-lid - 0.a58 


re ene ere ee en i ee eee eee a me m= me es ee ee ee 


we ee me 








ol) EXPERIMENTAL DETERMINATION OF TRANSIENT RESPONSE 

In order to determine the transient response of the motor, a step 
input was applied to the grids by inserting a double pole, double throw 
‘Switch across a variable voltage source. When the switch was thrown a small 
step input was applied. The output of the motor was recorded on the Miller 
oscillograph. A sample record is shown in Fig. 43. As can be seen from the 
sample record, it is very difficult to read exact values from such a picture. 
A comparison of the experimental and theoretical results is made in section 
5.3 and the experimentally determined transient response is plotted in 


fag. 10. 


{0. 








FIG.43 SAMPLE OF MILLER OSCILLOGRAPH RECORD OF TRANSIENT 
RESPONSE TO A STEP INPUT 





5.0 COMPARISON OF TREORETICAL AND EXPERIMENTAL RESULTS 
In this section the theoretical results obtained in section 3.0 


will be compared with those obtained in section }.0 by actual experiment. 


5.1 COMPARISON OF THEORETICAL AND EXPERIMENTAL IMPEDANCE CURVES 
In Fig. lh has been plotted the theoretical impedance curve for 

Motor B assuming a value of if =0O.1. In the same figure, the impedance 
curve obtained experimentally is plotted. Data for these curves are given 
in Table XIII. The agreement of the theoretical and experimental curves is 
fair except at low frequencies. The reason for this disagreement below 
resonance is not known but may be due to invalid data at the low frequencies 
or harmonic distortion since the data was taken from one coil only. The 
theoretical curve was calculated assuming Class A push-pull operation of the 
amplifier. Since the phase and amplitude curves fit so nicely in the next 
section, it is felt that the difficulty lies in the experimental techniques 


used in measuring the impedance curve. 
4 


5.2 COMPARISON OF TREORETICAL AND EXPERIMENTAL PHASE AND AMPLITUDE 

CURVES 

In Fig. 45 the theoretical amplitude characteristic for J =0.1 
and the actual experimental curve for Motor B are plotted for the purpose 
of comparison. The calculation data for these curves may be found in 
Tables XI and XIV. The correlation between the two curves is excellent. 

In Fig. 6 are plotted the phase characteristic curves correspond- 
ing to these amplitude curves. The agreement between theoretical and ex- 
perimental curves is remarkable as the curves would lie one on the other out 
to 60 cycles if we ignore the constant phase lag difference of 6 degrees of 
the experimental curve. This lag appears to be caused by the initial value 


of displacement of the armature (section lh.3) or by some delay in the phase 
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TABLZAAITII THAORETICAL (f=0.1) AND zXPERISENTsaL DATA FOR 
IMPhDANC. CURVzS OF AOTOR B 
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THSORETICAL EXPERIMENTAL ! 
| f zx10° Voltage(<7n) Current Voltage(coil) Lx19°. 
a ohms volts Be VOLES onws 
a 0 Beco) y 
/50 8.56 | .0% 2.67 40.0 15.00 
EO eS | One &.723 85.0 30.6 | 
P40 c8.1 |}. 8k 3.04 oes 34.6 | 
150 21.3 03% 3.08 joes 54.4 
160 $5.4 j| 031 Seow Poon 26.0 
yi79 = «40..5 074 2.74 ioe) 5605mm 
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oD 09.6 | 252 1.v6 36.0 18.4 | 
| £10 | a - 065 i240 Lies Te 
| Be ide. @ O74 ei Sono loa 
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THMORETICAL DATA FOR PHASE AND AsiPLITUDE FREQUE 
CiuARACTEHISTICS OF MOTOR b 
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measuring circuit or both. At 100 cycles per second the actual phase lag 
differs from the theoretical by only 12 degrees and the absolute phase lag 
is only 36 degrees. Beyond resonance there is more phase lag but this is 
accounted for by the fact that theoretically the system was assumed to be 
third order. For a third order system the phase lag should be asymptotic 
to the 270 degree axis. Actually, as shown in section 3.115, taking the 
elastance of the connection of the torque motor to the transfer valve into 
account the system is fifth order and therefore asymptotic theoretically 
to the 450 degree axis. 

It is interesting to check the value of effective mass of the 
system found in section .1. 


or M 


SP ligk =~ Kk 
nee M n° 


but i — Grp where &,is the frequency of maximum 
| fe Zi ~ response 
\ is damping ratio = 0.1 
The value of ©, from Fig. 45 is (27188) radians per second. Substituting into 
the formula for M, 
M. 168 = iets 1074 1b sec‘/in (20.8 grams). 


Gey 
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S.21 Survey of Linear and Non-linear Regions 





In performing the experimental tests for phase and amplitude 
characteristics on the torque motor, one must be careful to avoid voltage 
limiting, i.e. not allow the voltage on the tube to drop below approximately 
65 volts so that the tube goes into its non-linear region and thereby in- 
troduces non-linearities into the system. The tube characteristics are shown 
in Fig. 7. Expressed differently, one must exercise care not to drive the 
tube from pentode into triode operation. To ensure that in running the 


experimental tests data were secured that could be compared directly with 


8h. 
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FIG.47 TUBE CHARACTERISTICS OF THE 6AQ5 PENTODE 





that of the analytical analysis, an oscilloscope was placed across the tube 
aS a voltmeter. The voltage across the tube was kept at all times above 

65 volts. Also an initial value of amplitude of armature movement of 
t0.005 was used to avoid any non-linear effects due to the back electro- 
motive force, Kp- Avoiding these non-linearities accounts for the fine 


correlation between the theoretical and experimental data. 


5.3 COMPARISON OF THEORETICAL AND EXPERIMENTAL TRANSIENT RESPONSE 
CURVES 
In Fig. 8 is plotted the theoretical transient response to a 
unit step function. Also the experimental response to a step input with 
the amplifier connected and no 0.1 microfarad capacitors across the motor 
coils. The experimental curve was taken off the photographic record from 
the Miller oscillograph and was blown up proportionally to the scale of 
“the theoretical curve. The agreement is very good considering the dif- 


ficulty in reading values off the record. 








(q@ HOLOW) SHAYND ASNOdSHY LNAISNVUL TYINSMINSd Xa Any ieee oe HO NOSTYYd"00 a eed 


a | a ae ao 


pate 8) ig 


| 


Se 
iia 
i 


eter ee ee od re 








88. 


6.0 DISCUSSION OF EFFECT OF MOTOR STUDY ON CLOSED SERVO LOOP 
In this section the effect of the motor on the functioning of the 
closed servo loop as shown in Fig. 1, will be discussed. The transfer 


function of the motor itself from section 3.23 is: 


Ky/kR 
(s + 2950 S + 203) + j1190] [ (s + 203) — j1190] 





If we place the motor in series with an integrating device such as a hydraulic 


transfer valve and actuator with transfer function K/s, as shown in Fig. h9, 


MOTOR INTE GRATING DEVICE 











Ku/kR 
[S+ 29505 +203}j//FOK+203)~j1/ 90 


FIG. 49. BLOCK DIAGRAM OF MOTOR IN SERIES WITH INTEGRATING DEVICE 


we have for the open loop transfer function, 
K Ky/kR 
s [s+ 2950 |[ (s + 203) + j1190| | (3 + 203) - 51190] 
Substituting into this expression the values of the constants and multiply- 
ing out the quadratic term in the denominator, we may write the transfer 


function as 


K (9.64x107°) 
s(8 + 2950)(s¢ + 06s +1.h6x10°) 


In terms of time constants this can be expressed as 


1 where {= 
%,8(Ts +1)( (Gs)? + h06Gs +1) Kx2.23x107 
5 
2 2950 
1 


“3=T3T6 , Y=0.168 


iii sii ill alla 
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This open loop transfer function may be plotted by the use of non-dimensional 
gain phase diagrams which are included as Figs. 50 and 51. By making the 
plot, it was found that the locus crossed the minus 180 degree phase lag line 
at a value of @® equal to 115 radians per second. The table below shows the 
values of the individual transfer functions for a value of @ equal to 11)5 


radians per second. 





TRANSFER FUNCTION GAIN g 
db degrees 

EO -0.5 =20 

eee 10.0 -70 

wey o 

-180 


Since it is desired that the closed loop gain not exceed 3 decibels at minus 
180 degrees, the total open loop gain must not exceed minus l.5 decibels 

(see Fig. 52). It is evident then that the gain of the integrating neteenee 
@s, must equal minus 1.0 decibels. From Fig. 50, u equals 5.0 for a gain 


of minus 1.0 decibels. 


a G2 = 5.0 
2, 23x107 LK 
or k= 115 =102.5x101° volts/in-sec. 


2, 23x1071 x5 .0 


Therefore, K = 2,23x10719x102.5x10L° = 228/sec. 
Loop 


It can be seen that the motor allows extremely high values of gain 
to be used which is of great practical importance. In regard to this it is 
interesting to note the effect of increasing the damping term, J » in the 
quadratic term of the transfer function. For J =0.168 there resulted an 


open loop gain of 228; but, if j is increased to a value of 0.300, the open 
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loop gain, becomes equal to 456. This is closer to the values of 


Koop? 
gain actually realized in the closed loop shown in Fig. 1, since the 
hydraulic transfer valve adds a considerable amount of viscous damping to 


the system. 


6.1 EFFECT OF CURRENT AND VOLTAGE FEEDBACK 


<a 


Co, 
d R 
& 
g 
R 
= = Bt 


FIG. 53. SCHEMATIC SHOWING CURRENT FEEDBACK CIRCUIT 
To understand what happens if we put in a resistor, R? to 
ground forming a current feedback path, let us first write the voltage 


equation by Kirchhoff's law 


eg, = Ri + Ldi + Kpdx where R= ry, + R, (Section 3.113) 
at at 


The voltage input to the grid, e,, is the sum of Go) and Rois 


g? 
Cg — Cel + Roi or Cg) = eg — Roi 
Substituting this value for ©p1 into the initial equation, 


Meg = (R+ RoHS + Ldi + Kpdx 
dt dt 


The effect of current feedback is to increase the amount of resistance in 
the circuit by Ros which makes the tube look even more like a constant 
current source than before. 

Voltage feedback has just the opposite effect tending to reduce 


the resistance and making the tube appear more like a constant voltage source. 





9h. 


7.0 CONCLUSIONS 

This motor, operating on the polarized relay principle, provides 
high frequency response and by reducing the motorvalve time delay improves 
the servo stability, allows a higher servo loop gain, and results ina 
higher performance servo system. 

The high resonant frequency makes for good servo stability and 
keeps transient oscillations of the motor from appearing at the output of 
the servomechanism. 

Motors of this type will have many applications where small 
proportional displacements and large torques with low inertia are prime 


requisites. 
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8.0 APPENDIX 


8.1 MOMENT OF INERTIA OF ARMATURE CALCULATION 





FIG. 54. MOTOR ARMATURE DIVIDED INTO SIMPLE GEOMETRIC SHAPES 
The armature was divided up into simple geometric shapes as shown 


in Pig. Sli and listed below. 


Volume Mass i 2 iL 3 
Part Cu. In. Grams gm=cm gm-cm 
A (Rectangular 
parallelepiped) 2( 0.0927 248-95) 242526) “207 ailic) 
B (Circular 
cylinder) 0.1251 ¥6n7 --- 1.5) 
C (Rectangular 
parallelepiped) Oro) eng 0.26) 2. 26 
D (Circular 
cylinder) 2(0.0238) 6.13 eos are 
Total I, = 8.12 


Sample calculation for one body (rectangular parallelepiped A ) 
Volume = (0.6715 x 1.07 x 0.132) = 0.0927 cu. in. 


Mass =0.0927 cu. in. x 0.28 lbs/in? (density of steel)x Su gr/lb 
x 0.75 (staking factor) 


=8.95 grams 


I(Moment of inertia about axis h-h’)_ 8.95 (1.708°+ 0.3362) _. 2.26 gm-cm* 
12 
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Using the parallel axis theorem 
Iga, = T+md? =2.26 + (8.95 x 1.665) =17.16 gn-cmé 


This procedure was followed for each part of the armature and the result 


tabulated above. The moment of inertia of armature, I,, equals 48.12 gm—cm?. 


8,2 FLUX DENSITY IN AIR GAP DUE TO CURRENT IN ARMATURE COILS 
The reluctance of the air gap is 


& = = 0.15 — 0.0683 NI/maxwell 
He We 2a9al.09 


where g = length of gap, 0.15 cm. 
S = area of gap, 1.69 sq. cm. 
/4e = magnetic permeability, 
1.259 maxwell-cm/NI 


The flux in each air gap due to the current in the armature coils is 


[een NI, where NI, is the armature ampere- 
e e 
KR ' turns acting at each gap. 


The flux density is 


B. 0, — Mar_ NI, 
“S BRS  2x0.0683x1.69 


If we denote this flux density by the symbol Bga, then 


gauss 


Bga— NI, gauss where N = 7650 turns 
2x0.1152 Pga =flux density in gap due to 


armature, gauss. 

The measured values of the flux density were obtained by removing the 
permanent magnet from a motor and measuring the flux present in the gap for 
different values of current in the armature coils. The actual physical 
measurement of the flux density in the air gap presented quite a problem 
because of the small physical dimensions involved. A special test coil was 
wound by Mr. Loran A. Wenrich and imbedded in a piece of lucite. This search 

coil and the fluxmeter used were then calibrated together at the National 


Bureau of Standards and the combination was used in making all flux measure- 











x ‘ . 
— 2 
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rans ; 








8.35 TABLE XVI CONSTANTS OF MOTORS A AND B USED IN STUDY 





CONSTANT MOTOR A 








| MOTOR B 
“FR resistance (Retr), onus | 61800=S| = «61800 
| Re coil resistance, ohms 1800 | 1750 
| ry plate resistance, ohms 60000 60000 
L inductance, henries BS | Z0 
K, electromechanical coupling 
constant, lbs/aap 90.1 100 
k spring coefficient, lbs/in £O7 168 
| K}, back electromotive force, 
volt sec/in 9.0 9.6 
f/M daaping parameter /sec 300 £64 
M mass, lb sec“/in 1,.188x1074 1.188x107* | 
i critical damping ratio 0.113 0.100 | 
f viscous friction, lb sec/in 0.0354 0.0283 | 
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COMPLETE LIST OF SYMBOLS 


Constant 

Maximum displacement of armature from neutral 
Terminal point in phase shift network 

Special constant used in section 3.115 

Constant 

Flux density 

Flux density in air gap due to armature ampere-turns 
Flux density in any gap due to permanent magnet 
Terminal point in phase shift network 

Special constant used in section 3.115 
Capacitance 

Special constant used in section 3.115 

Terminal point in phase shift network 

Special constant used in section 3.115 

Voltage applied to armature coil 

Grid voltage 

Voltage input 

Voltage output 

Voltage 

Special constant used in section 3.115 

Viscous friction coefficient and also frequency 
(no ambiguity) 

Viscous friction of valve 

Force 

Force of valve 

Terminal point in phase shift network 

Special constant used in section 3.115 


Forces acting at air gaps 1, 2, 3, and 
Magnetic potential difference 
Gap length with armature centered 


Gap lengths of air gaps l, 2, 3, and 4 
Differential current 

Individual currents flowing in each of the two coils 
Current 

Moment of inertia of armature 

Current in armature coil 

Effective spring coefficient 

Anti-spring rate 

Actual spring rate of torsional centering spring 
Spring coefficient of valve 

opring coefficient of connection 

Gain of integrating device 

Back electromotive force 

Open loop gain of system 

Electromechanical coupling constant 

Step input voltage 

Length of armature suspension 

Natural logarithm 

Inductance of armature coil 

Total effective moving mass 





Rx 


Ry, R ,! 
R 


Re 


R1 Ro 434), 
8 


R 


t 


y, 


ERP H om 


AX 


~ 


Lor. 


Mass of armature 

Mass equivalent of armature 

Mass of valve piston and piston rod 
Number of turns on armature coil 
Ampere-turns of permanent magnet 
Ampere-turns of armature 

Pivot point of armature 

Magnetic permeance 


Permeances at air gaps 1, 2, 3, and 

Force of valve and rod acting on motor 

Distance from point of application of valve force 
to pivot point P 

Plate resistance 

Total resistance of plate and coil 

Resistance of armature coil 

Variable resistance 

Matched resistances in phase shift network 
Magnetic reluctance 

Reluctance of armature 

Reluctances at air gaps l, 2, 3, and 

Laplace operator 

Area of gap 

Period of oscillation 

Dimensioniess frequency variable 

Voltage generated by motion of armature in magnetic 
field 


Voltages 

Half suspension width 

Armature displacement at point of application of 
valve force 

Displacement of piston rod 

Length of armature from pivot point to center of 
air gap 

Capacitive reactance 

Armature displacement in air gap 

Two consecutive values of amplitude of same sign 
from oscillographic record of transient response 
Length of armature from center of air gap to point 
of application of valve force 

Impedance 

Impedance of armature coil 

Source impedance 


Logarithmic decrement 

Time interval 

Error signal 

Ratio of actual damping to critical damping 
External input signal 

Output signal 

Amplification factor 

Magnetic permeability 

Time constants 





SPEEPS 


noes 
£3, 24 
f1', Bo! 
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Phase angle 

Armature flux 
Leakage flux 

Flux at any gap 
permanent field flux 


Total flux at air gaps 1, 2, 3, and 


J3' £),'’ Loop fluxes used to calculate leakage flux 


Angular velocity 

Angular velocity at undamped natural resonant frequency 
Angular velocity at damped resonant frequency 

Angular velocity at peak resonance 

Ohms 


Bar over a quantity means the Laplace transform 

of that quantity 

Single dot over quantity indicates the first derivative 
of that quantity with respect to time 

Double dot over quantity indicates the second derivative 
of that quantity with respect to time 
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